Hamiltonians and Green's functions which 
interpolate between two and three dimensions 



Rainer Dick 

<N 

| Department of Physics and Engineering Physics, University of Saskatchewan, 

■ 116 Science Place, Saskatoon, SK S7N 5E2, Canada 

U 

<. 

■^J- . Abstract 

<N _ 

I propose to use Hamiltonians which contain two-dimensional and three-dimensional terms for 
the description of two-dimensional systems in physics. 

As a model system the evolution of three-dimensional wavefunctions in the presence of an 
I/*-) \ infinitely thin layer is studied. The model predicts distance laws for correlation functions which 

interpolate between two-dimensional and three-dimensional behavior. It also predicts that in 
certain cases transmission probabilities through thin layers should depend not only on the trans- 
| verse, but also on the longitudinal momentum of the infalling particles. 
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The model also yields a static potential which interpolates between the two-dimensional loga- 



rithmic potential at small distances and the three-dimensional 1/r-potential at large distances. 



1 Introduction 



X 

Two dimensions played a prominent role in the development of physics in the last twenty 
years. 

On the experimental side this was driven e.g. by the needs of very large scale integration, 
by applications of semiconducting layer structures, by exploitations of surface catalytic 
effects, and by the development of atomic-scale surface analysis and manufacturing tools 
like scanning tunneling microscopy and atomic force microscopy, to mention only a few 
developments in this 

On the theoretical side interest in two-dimensional field theories was largely driven by 
string theory [1] , where fundamental excitations are described by covariant two-dimensional 
field theories, and by the realization that in two-dimensional critical systems with a ro- 
tational symmetry scaling symmetry may be elevated to full conformal invariance [2]. 

The purpose of the present paper is to point out that recent developments in the 
mathematical formulation of brane world models may also inspire new developments in 
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the physics of low-dimensional systems, and help us to acquire a better understanding of 
the transition between three-dimensional and two-dimensional behavior in these systems. 

The present work was specifically motivated by the brane world model of Dvali, Gabada- 
dze and Porrati, who recently proposed and analyzed a model which combined gravity on a 
(3 + l)-dimensional manifold (a "3-brane") with gravity in an ambient (4 + l)-dimensional 
bulk [3]. They observed that the combination of gravity in different dimensions yields 
a gravitational potential which interpolates continuously between the three-dimensional 
—1/r potential at small distances and the four-dimensional — 1/r 2 potential at large dis- 
tances, with a transition scale Idgp — rn\/m\ set by the ratios of the reduced Planck 
masses on the brane and in the bulk. Since we are not concerned with gravity in the 
present paper I will not explicitly write down the model in terms of intrinsic and extrinsic 
curvature terms, see [4], but instead refer to the related model of Dvali, Gabadadze and 
Shifman [5], which combines a Maxwell term in Minkowski space (coordinates x M = {£, r}) 
with a Maxwell term in an ambient (4+l)-dimensional bulk (coordinates x = {t, r, x 1 }): 



S = -^jdtjd^F, v F^ 



x ± =0 



1 



dt /cfr I dx ± FmnF 



The resulting Coulomb potential on the (3 + l)-dimensional Minkowski spaced is [5] 
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with the sine and cosine integrals 
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The ratio of gauge couplings defines a length scale 
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and A interpolates between a three-dimensional distance law at short distances and a 
four- dimensional distance law at large distances: 
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Action principles of the kind (1) were denoted as dimensionally hybrid action principles 
in [6]. 



Sec Rcf. [6] for a discussion of dynamical potentials on the brane and in the bulk. 
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Of course, this does not simply carry over to low-dimensional systems in condensed 
matter or statistical physics: Dimensionally hybrid action principles would not be a suit- 
able tool for model building in theoretical investigations of these systems because time 
derivatives generically will appear only as bulk terms in the Lagrangian of a particle 
interacting with a low- dimensional structure. 

Therefore the main proposal of the present work is to use dimensionally hybrid Hamil- 
tonians in model building for low- dimensional systems: A combination of two-dimensional 
and three-dimensional terms and the ensuing interpolating correlation functions may help 
to narrow the gap between the powerful methods of two-dimensional field theory and re- 
alistic thin layers or surface structures in physics and technology 

In the sequel I will use this idea to discuss non-relativistic particles interacting with a 
thin layer. The system becomes a dimensionally hybrid system with a specifically two- 
dimensional component through the assumption that particles in the layer have a kinetic 
energy different from particles outside of the layer, e.g. as a consequence of mass renor- 
malization M — > m due to the interaction of the particles with the components of the 
layer. This yields some straightforward but interesting results. 

In the next section I will show in a simple model that transmission probabilities through 
thin layers in this class of models depend also on the momentum parallel to the layer if 
\x = m(L) I X|i_+Q remains finite. The discussion of Green's functions in these models will 
be the subject of Sees. 3 and 4. Sec. 3 contains in particular a potential which interpolates 
between two-dimensional and three-dimensional distance laws. 



2 Dimensionally hybrid Hamiltonians and Green's functions 

To investigate implications of dimensionally hybrid Hamiltonians for the description of 
the interaction of particles with a thin layer we assume in the present section that the 
layer is planar and homogeneous and therefore generates a potential U(z), where z is 
transverse to the layer. 

In realistic two-dimensional systems particles are not strictly bound to the layer, and 
the effective particle mass in the layer may be changed due to interactions. This motivates 
the following Hamiltonian for particles of mass M in the presence of the layer 



h 2 r 

H=— d 2 xV^ + • W> 
2u J 




(4) 



z=0 



Here and in the sequel all vectors are 2-dimensional vectors in the layer. 

The assumption behind the Hamiltonian (4) is that the same field ip may describe e.g. 
free electrons in the bulk and large polaronsj^] or other collective excitations involving 
conduction electrons in the layer. The parameter /i has dimensions of mass per length, 
and in a limiting procedure from layers of finite thickness L would correspond to 



See e.g. [7,8,9] for introductions to polarons in solids. 
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(5) 



where m(L) would be the mass of the modes in the layer. 

The corresponding equation of motion for stationary single-(quasi)particle wavefunc- 
tions is 

ity(x, z) = -6(z)^-Ai>(x, 0) - ^ (A + d 2 ) ^(x, z) + U(zMx, z). (6) 
The Fourier ansatz 

■0(x, z) = — [(fkipik, z) exp(ik • x) (7) 

2,71 J 



yields the separated equation 

E ~ w) m z) = ~S^ (k ' z) + u{zMK z) + 6{z) ^ {K 0) - (8) 



Every solvable model of one-dimensional quantum mechanics gives a solution to this 
class of layer models, with the kinetic term of the layer modes only generating a cusp 
proportional to (M/fx)k 2 in ln^(k, z). 

Obviously, the large longitudinal momentum modes are strongly affected by the exis- 
tence of layer modes, but we will see in a moment that the two-dimensional kinetic term 
can also have a strong impact on modes with small longitudinal momentum. 



Fig. 1: Contribution of a virtual planar mode to a particle penetrating a thin homogeneous 
layer. 

It is a trivial exercise to adapt solutions of one-dimensional quantum mechanics to the 
cusp imposed by the layer modes, but it may be worthwhile to write down the modifi- 
cations of the transmission coefficient due to the layer modes when the layer potential 
represents a work function: U(z) = —w5(z). The transmission coefficient for an infalling 
particle of momentum {k, k±} is 



M 2 ( w 



(9) 
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i.e. the layer modes increase the transmission probability for low longitudinal momentum 
modes < h 2 k 2 < Afiw, and decrease the transmission probabilities for the modes of 
higher longitudinal momentum. 

This model also predicts a resonance in transmission for a certain value h 2 k 2 = 2fiw 
of the longitudinal momentum. This is as a genuine consequence of the two-dimensional 
kinetic term in (4) and may be the simplest way to test the viability of the idea of 
dimensionally hybrid Hamiltonians in low- dimensional systems. 

Obviously, the prediction of dependence of transmission probabilities on longitudinal 
momenta requires finiteness of the parameter \i (5), i.e. a derivation of the "phenomeno- 
logical" Hamiltonian (4) from a limiting procedure of purely three-dimensional models 
will require a thorough study of finite size effects on mass renormalization in solids. 



3 Correlations on a layer 

A model similar to (4) allows for a neat discussion of the impact of combinations of 
kinetic terms from different dimensions on the "two-dimensional" correlation functions on 
the layer. 

For this we assume that the layer is not necessarily homogeneous, but generates a 
strongly localized potential 

V(x,z) =u(x)5(z) 

along the layer: 

H= — Jd 2 *Jdz (W + • + d z ^ + ■ d z ip) (10) 

r ( h 2 

+ J d 2 x I — • Vip + ip + ui> 

The generating functional for correlation functions on the layer 

Z[j,J + ] = j#^ + exp(-/3tf[^+] - |rf 2 x[^+(x,0)j(x) +j + (x)^(x,0)]) (11) 

= ex P (-/3/rf 2 x^n(x)^_) Z \j,j+} 

with 

Z \jJ + ] oc exp^ / rf 2 x / dVj + (x)G(x - x', 0)j(x')) • (12) 

Since (10) is a free theory from a field theory point of view, the "two-dimensional" corre- 
lations in it can be calculated from tree-level diagrams, which involve only the restriction 
of the free Green's function G(x — x', z) to the layer and insertions of the layer potential. 

4 As usual S/5j acts from the right if i/j is fermionic. 
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The free Green's function used in (12) satisfies 

(A + df) G(x, z) + 2£S(z)AG(x, 0) = -S(x)6(z) (13) 
where 

M 

21=—. (14) 
fi 

The ansatz 

G(x, z) = — — J d 2 k J dk ± G(k, jfej.) exp[i(k ■ x + k ± z)] (15) 
yields 

(k 2 + k 2 ± ) G(k, fc ± ) + ^A; 2 /" d/^ G(k, fc^) = 1. (16) 
This determines the k±- dependence of the propagator 



G ^) = 4%- (17) 



With 
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we find the Green's function 

1 



G(k, fcj 



(l + JW)(Jfe 2 + ifei) v 7 

or 

G < k ' z > = S(iW expHcW) ' (19) 

The solution^ 

G( x„) ■^7 a 7 ^^-M)i .^; a s^p w>r) (20) 



of (13) can be thought of as the electrostatic potential of a unit charge on the layer, if the 
fields which are continuous across the layer make a special contribution to the Hamiltonian 



I follow the conventions of [10] for Bessel and Struve functions. 



of the electromagnetic fieldQ: 

H[F]=i Jd 2 * (E 2 + 5i) + ^y>x Jdz (E 2 + E\ + B 2 + B 2 ± ) 
e.g. as a consequence of a non- vanishing limit of 



(21) 



21 = lim(e r L). 



Here e r is the relative permittivity of the layer and L its transverse extension. 

The perturbation series (11,12) requires the Green's function on the layer, which can 
be expressed as a linear combination of a Struve function and a Bessel function of the 
second kind: 



G(x, z) 
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This interpolates between two-dimensional and three-dimensional distance laws 
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$(x) along with the limiting cases is plotted in Fig. 2. 



Fig. 2: The solid line is the Green's 
function (22) on the layer as a function 
of x = r/£, in units of t~ x . The up- 
per dashed line is the three-dimensional 
l/A.Tir potential, and the lower dashed 
line is the two-dimensional logarithmic 
potential. 
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6 In the spirit of the philosophy advocated in this paper the electromagnetic fields in (21) would 
simultaneously describe photons in the bulk and polaritons in the layer. 
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Polaron masses in semiconductors are of the order of effective band masses. Due to 
strong band curvature the ratios of polaron masses to the mass of free electrons can be of 
order m/M ~ 10~ 2 , so one may hope that in thin semiconducting layers of thickness L the 
parameter i is of order t ~ 10 2 L. Two-dimensional distance laws for correlation functions 
might then be realized up to distances of order 10L, and intermittent behavior for distances 
between 10L and 200L. For a possible realization of (22) as an electromagnetic potential 
in thin layers, semiconducting compounds involving Pb might be the best bet due to their 
high relative permittivities of order e r ~ 10 2 — 10 3 [11]. The corresponding field strength 
on the layer is 



- <9 r $(x) 
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Hi 



-Y 1 



(23) 



Fig. 3: The solid line is the field 
strength per charge (23) on the layer as 
a function of x = r/£, in units of l~ 2 . 
The dashed line approaching the solid 
line for x < 1 is the two-dimensional 
l/Anir field, and the other dashed line 
is the three-dimensional l/4nr 2 field. 
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This reduction of the force between charges in a thin dielectric layer with finite 2£ = 
lim i _ +0 (e r L) can pictorially be understood as a consequence of the fact that field lines 
are refracted away from the layer when they leave the layer. This reduces the field lines 
e.g. between two opposite charges at short distances, since the field lines cannot re-enter 
the layer on short scales, whereas for large separation of the two charges along the layer 
re-entry renders the refraction effect negligible. 



4 The Green's function for scattering from two-dimensional potentials on 
the layer 



The stationary wave equation from (10) is 

2fl 



z) = 5(z) [ ~A + w(x) J V(x, 0) - ^ (A + d 2 z ) ^(x, z). 



(24) 



S 



We have already noticed that this can be solved exactly for w(x)=constant. 

In discussing scattering of bulk particles from the layer in the model (10) we could 
proceed using the ordinary three-dimensional Green's function for scattering of waves of 
energy E = h 2 K 2 /2M and treat the full two-dimensional contribution to Eq. (24) as a 
perturbation. However, here I rather would like to treat only the layer potential -u(x) as 
a perturbation. This has the virtue of reducing the perturbation for large longitudinal 
momenta. 

The relevant unperturbed wave is then 



^o(x, z) = — l=g exp(iK|| • x) 

V 2"7T 



T-f? tPT~d ex P( — l K±z) 



(25) 



+ ^ Z) K ± +\ K fi eMiK±z) 

where again the definition (14) was used. 

The relevant Green's function Gk for propagation of bulk plane waves of energy 

2M 2M V II ± J 
has to satisfy 

(A + d 2 z + K 2 ) G K (x, z) + 2£5(z)AG K (x, 0) = -6(x)6(z), (26) 
and the solution proceeds similarly to the solution of (13). The Fourier ansatz (15) yields 

(k 2 + k 2 ± - K 2 ) G K (k, k ± ) + Uz 2 j dk' ± G K (k, k' ± ) = 1, (27) 
which determines the /c^-dependence of the propagator 

°*™ = X> + f ki-K> - (28) 

At this stage the possibility of poles complicates the calculation slightly: In evaluating 
the integral in (27) with (28) for k < K we have to make a judicious choice on how 
to shift the poles or the integration path at k± = ±V K 2 — k 2 , corresponding to correct 
physical boundary conditions on GkQs.,z). The correct choice turns out to be k± = 
±(\/K 2 — k 2 + ie) since G K (k,z) is supposed to describe outgoing scattered waves from 
the layer if k < K, i.e. we have 

g '^ = * + *i-V fe (29) 
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With 



1 °° 
- / dk' ± 

71 J 



Q(k 2 -K 2 ) + .Q(K 2 -k 2 ) 



k 2 + kf- K 2 - ie Vk 2 - K 2 VK 2 - k 2 

—00 

we find the Green's function at large longitudinal wavelength k < K 

\JK 2 — k 2 

G ^ Z ^ 2k*(-2WK-<-^ XP ( lVW ^^' <31) 
while the short longitudinal wavelength part is 

G K (k, k ± ) = ^= _ R2 + k2 Q~ 2 + k 2 ± _ R 2_ ie) ' ( 32 ) 

G K (k, z) = 1 exp(-Vk 2 - K 2 \z\) . (33) 

y 1 2Vk 2 - K 2 + 2kH V 1 y ! 

Of course, the Green's function again reduces to the usual three-dimensional result for 
£ — > 0, i.e. if the modes in the layer become so heavy that they decouple. 
With an incoming plane wave the integral equation following from (24,26) is 

2M r 

V>(x, z) = V>o(x, z )--fpJ rfV G ^ - x '' ^(x')V(x', 0). (34) 

In a Born approximation this yields with (25) 

2 M K f 

V>(x, z) = ^o(x, 2) 7 — — i 2 / ci 2 k exp(ik • x)G^(k, z)u(k - Ky), (35) 

V27T ft 2 A -L + lK \\Z J 

where the normalization of the Fourier transformed layer potential is 

u(q) = J d 2 x. exp(— iq • x)-u(x). 

Eq. (35) together with (33) implies that no particular scattering wave is generated by 
the short wavelength components at |q| < K\\ + \J K 2 + K\ of the layer potential, which is 
simply a statement of the limited resolving power of the external wave. However, beyond 
that the ^-dependence of (31) and (35) implies that the two-dimensional kinetic term in 
(10) reduces potential scattering at large K\\. 

5 Summary 

Two-dimensional field theory is a very appealing subject with many powerful results. 
One reason for this is because every tensor and spinor field on a 2-manifold decomposes 
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into covariant primary fields which provide 1-dimensional representations of the corre- 
sponding symmetry groups[3]. 

However, the assumption of strict two- dimensionality seems too restrictive when it 
comes to comparisons with actual layers or surface structures in physics. Conservative 
theoretical approaches to low- dimensional structures in physics and technology therefore 
rely on genuine three-dimensional Hamiltonians and only restrict the locations and mo- 
menta of particles to a surface or a layer (see e.g. Sec. 9.2 in [8]). In these approaches 
two-dimensionality is only taken into account at a kinematical level, at the expense of 
sacrificing the powerful methods and results of two-dimensional field theory. 

On the other hand, recent results in brane theory taught us that straightforward com- 
binations of four- dimensional terms and five-dimensional terms in action principles yield 
interpolating Green's functions on the brane, and it is apparent from the functional inte- 
gral representation that this property must also hold for higher order correlation functions 
on the brane. 

As mentioned above this cannot carry over directly to low-dimensional systems in 
condensed matter and statistical physics, but it initiated the present proposal to use 
dimensionally hybrid Hamiltonians H = Ihi + H3 for theoretical investigations of low- 
dimensional structures in physics. 

Such an approach has the prospect to provide more realistic results than strictly two- 
dimensional field theory, while at the same time utilizing the power of two-dimensional 
field theory for the determination e.g. of equilibrium correlation functions in the limiting 
cases £ — > 00 or k£ 3> 1. 

Further virtues of this approach are predictions on the transition behavior between two- 
dimensional and three-dimensional distance laws in layers and surface structures, and a 
better understanding on how two-dimensional structures might be approached in the more 
conservative purely three-dimensional framework, through studies of finite size effects on 
effective masses and permittivities in three-dimensional models. 

To illustrate the use and some straightforward consequences of dimensionally hybrid 
Hamiltonians a homogeneous layer and layers with strongly localized potentials were 
studied. In these settings the field ip describes simultaneously bulk particles of mass M 
and excitations of mass fi per transverse length in the layers. These models might serve 
as an approximation to thin layers of semiconductors or polar solids, where the planar 
modes would correspond to large polarons. The model predicts a strong dependence of 
transmission probabilities on longitudinal momentum. 

The Green's function in the model also yields a static potential which interpolates 
between the logarithmic two-dimensional distance law at distances <C i and the three- 
dimensional Coulomb law at distances ^> i. 

Acknowledgement: This work was supported by NSERC. 



7 This is true beyond the realm of conformal transformations if the Beltrami parameters on the 
2-manifold are used to decompose tensors and spinors into covariant primary fields, see Sees. 1 
and 2 in [12] for tensors and [13] for spinors. 
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